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Abstract 

In this paper we study the production of a heavy charged Higgs boson in association with heavy 
quarks at the LHC, in a type-H two-Higgs-doublet model. We present for the first time fully- 
differential results obtained in the four-flavour scheme at NLO accuracy, both at fixed order and 
including the matching with parton showers. Relevant differential distributions are studied for two 
values of the charged boson mass and a thorough comparison is performed between predictions 
obtained in the four- and five-flavour schemes. We show that the agreement between the two 
schemes is improved by NLO(-I-PS) corrections for observables inclusive in the degrees of freedom 
of bottom quarks. We argue that the four-flavour scheme leads to more reliable predictions, thanks 
to its accurate description of the bottom-quark kinematics and its small dependence on the Monte 
Carlos, which in turn is rather large in the five-flavour scheme. A detailed set of recommendations 
for the simulation of this process in experimental analyses at the LHC is provided. 
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1 Introduction 

Charged Higgs bosons appear in several extensions of the scalar sector of the Standard Model (SM). 
In particular, as the SM does not include any elementary charged scalar particle, the observation 
of a charged Higgs boson would necessarily point to the presence of new physics. 

In this paper we focus on one of the simplest extensions of the SM featuring a charged scalar, 
namely the two-Higgs-doublet model (2HDM). Within this class of models, the existence of five 
physical Higgs bosons, including two (mass-degenerate) charged particles , is foreseen. Imposing 
flavour conservation, there are four different ways to couple the SM fermions to the two Higgs 
doublets. Each of these four ways of assigning the couplings gives rise to a different phenomenology 
for the charged Higgs boson. Here we consider the so-called type-H 2HDM, in which one doublet 
couples to up-type quarks and the other to the down-type quarks and the charged leptons. The 
Minimal Supersymmetric Standard Model (MSSM), up to SUSY corrections, belongs to the type-H 
2HDM category. 

Light charged Higgs scenarios are classified by Higgs boson masses smaller than the mass of the 
top quark (typically rnH± ^ 160 GeV), where the top decay to a charged scalar and a bottom quark 
is allowed. One refers to heavy charged Higgs bosons, on the other hand, for masses larger than the 
top-quark mass (typically > 180 GeV). In this case, the dominant production channel 

at the Large Hadron Gollider (LHC) is the associated production with a top quark/antiquark and 
a (possibly low transverse momentum) bottom antiquark/quark. In the intermediate-mass range 
(160 ;< 180 GeV) width effects become important and the full amplitude forpp —>■ H^W^bb 

(including non-resonant contributions) must be taken into account for a proper description. The 
intermediate-mass range has not been studied so far at the LHC Run 1. 

Searches at LEP have set a limit mH± > 80 GeV on the mass of a charged Higgs boson for 
a type-H 2HDM [^. The Tevatron limits have been superseded by results of the LHC 

experiments. Recent ATLAS results 0 for a type-H 2HDM based on 19.5 fb“^ of pp collisions at 
8 TeV exclude BR(t — bH'^) ■ BR(iJ^— )■ larger than (0.23 — 1.3)% in the low-mass region 

(80 GeV < mH± < 160 GeV). For the first time limits on tH'^ production times BR(iL'*' —>■ 
are provided for 180 GeV < < 1000 GeV. Rates above 0.76 pb at low masses and 4.5 fb at 

large masses are excluded. A preliminary note by GMS [^, based on the analysis of 19.7 fb“^ oipp 
collisions at 8 TeV in the same search channel, sets an upper limit on BR(t —)■ bH'^) ■ BR(i7+ —>■ 
t^Vt) between (0.16— 1.2)% in the low-mass region (80GeV < mH± < 160GeV). This limit 
supersedes the one published in Ref. [^. In the high-mass region (180GeV < < 600GeV), 

limits on production times BR(i/^ —>■ t^Vt) are set between 0.38 pb and 26 fb. Finally, CMS 
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Figure 1: Leading-order diagrams for heavy charged Higgs and top associated production, in the 
4FS (left and centre) and 5FS (right). 


has also published a preliminary note [H on a direct search for a heavy charged Higgs which decays 
in both the —>■ tb and the —)■ 7 ^ Vt channels. 

In this work, we consider heavy charged Higgs boson production at hadron colliders and leave 
the intermediate-mass range to future studies [^. In particular, we focus only on the production of 
a negatively charged scalar since the results are identical for a positively charged scalar. As for any 
process involving bottom quarks at the matrix-element level, two viable schemes exist to compute 
the production cross section of a heavy charged Higgs boson. These are usually dubbed as four- 
and five-flavour schemes. In the four-flavour scheme (4FS) the bottom quark mass is considered as 
a hard scale of the process. Therefore, bottom quarks do not contribute to the proton wavefunction 
and can only be generated as massive final states at the level of the short-distance cross section, 
entailing that 6-tagged observables receive contributions starting at leading order (LO). In practice, 
the theory which is used in such a calculation is an effective theory with four light quarks, where 
the massive bottom quark is decoupled and enters neither the renormalisation group equation for 
the running of the strong coupling constant nor the evolution of the parton distribution functions 
(PDFs). The LO partonic processes in the case at hand are 

gg —)• H~bt and qq —>■ H~bt. (1) 

Next-to-leading order (NLO) calculations for the total cross sections in this scheme have been 
presented in Refs. [^[lo]. 

Conversely, in five-flavour scheme (5FS), the bottom quark mass is considered to be much 
smaller than the hard scales involved in the process. The simplest definition of the 5FS—that 
suits particularly well perturbative computations—is to strictly set mb = 0 in the short-distance 
cross section. Consequently, bottom quarks are treated on the same footing as all other massless 
partons. The only difference is the presence of a threshold in the bottom-quark PDF and the initial 
condition of the bottom quark evolution being of perturbative nature. The use of 6-PDFs comes 
along with the approximation that, at leading order, the massless b quark has a small transverse 
momentum. In this scheme, the leading logarithms associated to the initial state collinear splitting 
are resummed to all orders in perturbation theory by the Dokshitzer-Gribov-Lipatov-Altarelli- 
Parisi (DGLAP) evolution of the bottom densities. The LO partonic process is given by 

gb^H-t. ( 2 ) 


Next-to-leading order predictions for heavy charged Higgs boson production in the 5FS 
including the matching to parton-shower Monte Carlos, were studied in Refs, 
corrections 


II 17 


possibly 
Electroweak 


18,19 and soft gluon resummation effects 15 20 21 have also been included in recent 


works. 

The leading-order diagrams in the 4FS and 5FS are displayed in Fig.j^ The comparison between 
the two schemes at the level of total cross section has been performed by several groups, see e.g. 
Ref. and references therein. In a more recent study 22 a thorough combination of all sources 


of theoretical uncertainties is performed, state-of-the-art PDF sets are used, the new scale-setting 
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23,24 is adopted and a Santander-matched prediction 25 


procedure proposed in Refs, 
four- and five-flavour scheme computations is provided. Furthermore, a comparison of five- 


of the 
and 


six-flavour schemes, assessing the effect of the inclusion of top PDFs, has been recently performed 
at the level of total cross sections 26 . In summary, these studies establish a generally good 


agreement of total cross section predictions in different flavour schemes, once a judicious scale 
choice is made. 

Studies at differential level of heavy charged Higgs production in four- and five-flavour schemes 
are very limited, if not absent. NLO plus Parton Shower (PS) accurate predictions are available 
only in the 5FS for both MC@NLO 16 and POWHEG [T^, and their comparison to 4FS results 


at the level of differential distributions has never been performed to date. Such a differential study 
Is certainly more involved than the inclusive one, mainly for two reasons. First of all, there is a very 
large number of possibly relevant observables at the differential level. Second, differences between 
the two schemes are generally larger than In the Inclusive case. As a consequence, an even more 
careful assessment of the related uncertainties is necessary for distributions. Indeed, observables 
Inclusive in the degrees of freedom of the bottom quarks, such as total cross sections, turn out 
to be quite similar In four- and five-flavour schemes. On the one hand power-suppressed terms 
are small for this class of observables. On the other hand collinear logarithms are phase-space 
suppressed 23 and therefore moderate, unless a very heavy charged Higgs Is produced. 

As far as exclusive observables are concerned. It is vital to investigate 4FS and 5FS predictions 
and assess their differences on a case-by-case basis in an unbiased manner, which is one of the 
main goals of this paper. The generic rationale is the following: if power-suppressed terms are 
relevant, the 4FS provides a more reliable prediction, while if collinear logarithms are large, the 
reorganisation of the perturbative series in the 5FS improves the stability of the perturbative 
expansion. 

The comparison of the two schemes is further complicated by the fact that observables related 
to light and 6-fiavoured jets are associated with different perturbative orders: while In charged 
Higgs production one final-state b quark (not accounting for top decays) Is already present at LO 
in the 4FS, a final-state b quark enters the 5FS computation only at NLO. Furthermore, tagging b 
quarks In a 5FS fixed-order computation necessarily leads to unphysical results, since they can be 
considered only as part of jets to retain infra-red cancellations. This issue clearly does not affect 
the 4FS, in which observables related to the b quarks are regulated by the b mass. Nevertheless, 
when considering b quarks at very large transverse momenta, 6-jet observables have to be preferred 
anyway, otherwise large logarithms \og(jpT/mb) would spoil the perturbative convergence. 

Matching to parton showers improves most notably the 5FS predictions, since the parton back¬ 
ward evolution of the initial-state bottom quarks turns them into massive final states and further 
generates 6-flavoured hadrons which renders 6-tagging realistic. One should bear in mind, however, 
that the details and the actual Implementation of such backward evolution are highly non-trlvial— 
being based on DGLAP evolution equations, which are only LL accurate—and therefore turn out 
to be widely Monte Garlo dependent (see Sect. 3.3 of Ref. 27 , and Sect. |3.3| below). Moreover, the 
necessity of reshuffling the massless into massive bottom quarks may have significant effects on the 
kinematics of final-state B hadrons. In the 4FS, the shower primarily Improves the description of 
Sudakov-suppressed small-pj’ radiation, by resumming leading collinear logarithms to all orders. 
In both schemes the PS Introduces additional power-suppressed terms in the soft region. 

In this paper, we present for the first time NLO-I-PS accurate 4FS predictions and thoroughly 
compare them to the 5FS ones. Our primary aim is to acquire a detailed understanding and 
assess which scheme is better suited for the simulation of the charged Higgs production signal, 
whose search represents a central part of the physics program in (and beyond) Run H at the LHG. 
As will be shown, our analysis confirms and extends the conclusions that have been drawn in 
previous studies of differential distributions In the four- and five-flavour schemes, in the context 
of Higgs production in association with bottom quarks 
production 


28 


29 . Similar analyses were performed for single top production 30 31 


and Higgs and single-top associated 

Moreover, 


given that the mass of the charged Higgs boson is still allowed to take almost arbitrarily large 
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values, this process is perfectly suited for a case-study of different flavour schemes in QCD in a 
broad range of scales. 

The paper is structured as follows: in Sect. an outline of the calculation is given and details 
of the implementation of the 4FS calculation at NLO-I-PS accuracy are provided. Results are 
displayed in Sect. in which the features of the 4FS calculation are described, including the 
size of the interference term that appears in this scheme, and the 4FS and 5FS distributions are 
compared. We conclude in Sect. with our final recommendations for experimental analyses. 


2 Outline of the calculation 


Our computation takes advantage of a full chain of automatic tools developed to study the phe¬ 
nomenology of new physics models at NLO QCD accuracy in the MadGraph5_aMC@NLO 
framework. 
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2.1 Framework 

Besides the usual tree-level Feynman rules, some extra ingredients have to be provided to the 
matrix-element generator in order to obtain the code for the simulation of a new physics pro¬ 
cess at NLO. These extra ingredients are the Ultra-Violet (UV) renormalisation counter-terms 
and a sub-class of the rational terms that enter the reduction of the virtual matrix elements, the 


so-called i ?2 terms 33 . UV and R 2 counter-terms can be computed starting from the model La- 
grangian for any renormalisable theory via the NLOCT package 34 , based on FeynRules 35 
and FeynArts [^. The tree-level and NLO Feynman rules are exported in the Universal Feyn- 
rules Output (UFO) format, as a Python module which can be loaded by matrix-element 
generators, such as MadGraph5_aMC@NLO. Finally, the UFO information is translated into 


helicity routines (38| by ALOHA 39 


We remind the reader that MadGraph5_aMC@NLO is a meta-code that automatically gen¬ 
erates the code for simulating any process at NLO(-fPS) accuracy. MadGraph5_aMC@NLO 


adopts the FKS method 40,41 for the subtraction of the singularities of the real-emission matrix 
elements, as automated in the MadFKS module [4^; virtual matrix elements are computed via the 


MadLoop module 43 , which relies on the OPP |44| and Tensor Integral Reduction (TIR) 46 
methods, as implemented in CutTools 47 and IREGI 48 , supplemented by an in-house im¬ 


plementation of the OpenLoops procedure 49 ; finally, the generation of hard events and their 


matching to parton-shower simulations is performed a la MC@NLO [^. The matching to Her- 
wig6 51 , Pythia6 (ordered in virtuality or in transverse momentum, the latter only for 


processes with no light partons in the final state), HerwigH — h 53 and Pythia8 54 is available. 
As a consequence, the only inputs needed are the implementation of the model in FeynRules and 
the inclusion of the running of the bottom quark mass, which follows from our renormalisation- 
scheme choice. More details are given in the next section. Both the 4FS and the 5FS computations 
have been performed in the framework specified above, with the additional advantage of ensuring 
the consistency of all settings and input parameters. 

To guarantee the validity of our analysis, we have performed a detailed comparison of the in¬ 
clusive cross sections obtained with MadGraph5_aMC@NLO against previous results published 
in literature. The calculation in the 4FS has been compared to the private code of Ref. 10 at the 


level of total rates, and a good agreement within the numerical uncertainty of the reference code 
has been found. Furthermore, with the bottom quark Yukawa renormalised on-shell and with its 
mass set the value of the top pole mass, we reproduce the tiH total cross section in the SM at the 
few per-mille level. 
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2.2 Implementation 


We have used the implementation of the generic 2HDM in FeynRules detailed in Ref. 34 . This 


model has been converted into a type-II 2HDM by adding jS as an external parameter and by 
restricting the Yukawa couplings accordingly. If top and bottom quarks are assumed to be the 
only massive fermions, the only non-zero entries of the Yukawa coupling matrices to the doublet 
without vacuum expectation value in the Higgs basis for the type-II 2HDM are given by 


Gg 3 = —cot j3 and Gg g = tan j3, 


(3) 


where are the Yukawa masses of the top and bottom quark. The parameter tan/3 = V’ljvx 
is the ratio of the vacuum expectation values Vi and V 2 of the two Higgs doublets, such that 
iP' = -\- V2 = ^ is the SM Higgs vacuum expectation value, where Gp is the Fermi 

constant. With those restrictions, the H~tb vertex is given by 


^tbH- — 


ytPR 


tan/3 


ybPLtan/3 


(4) 


where Pr/l = (1 ± 75)72 are the chirality projectors and yt/b = V^—^ are the corresponding 
SM Yukawa couplings. We strictly separate the Yukawa masses that are used in the computation 
of the couplings between the fermions and the scalars from the kinematic masses that are used 
everywhere else and set to the on-shell massQ This distinction allows us to keep a non-vanishing 


bottom Yukawa in the five-flavour scheme as the leading term in the small mb expansion IRI 14 


Furthermore, it allows us to choose different renormalisation schemes for the bottom quark mass 
in the matrix element and in the Yukawa coupling. 

The model i?2 and UV vertices required for NLO computations in MadGraph5_aMC@NLO 
have been computed using NLOCT 34 . The masses and the wave functions are renormalised 
in the on-shell scheme to avoid the computation of loops on external legs. The strong coupling 
constant is renormalised in the MS scheme with the contribution of massive quarks subtracted 
from the gluon self-energy at zero-momentum transfer. Therefore, only the massless modes affect 
the running of a^- The renormalisation of the masses in principle fixes the renormalisation of the 
top and bottom Yukawa since 

:6mt/b 


Syt/b = V2- 


with 


Smt/b = - 


9s 


) '^tjb 




(5) 


( 6 ) 


127r2 'ye 

in the on-shell scheme. This is the default renormalisation used in NLOCT and it would ensure 
that nothing but the strong coupling constant depends on the renormalisation scale. The top mass 
and Yukawa are always renormalised in this way throughout this paper. Therefore its Yukawa mass 
is set equal to the pole mass. On the contrary, the bottom quark Yukawa has been renormalised 
in the MS scheme, i. e. 


Syb = - 


V2 


1 V 

9 8^1 


(7) 


V 47r^e 

This scheme choice has the advantage of resumming potentially large logarithms \og{yLp/mb) (with 
yip ^ mH±) to all orders. The bottom Yukawa mass is set to the value of the running MS mass at 
the renormalisation scale. Besides the modifications at the level of the UFO model, also the code 
written by MadGraph5_aMC@NLO had to be changed in order to account for the additional 
scale dependence introduced by the 5-quark Yukawa, in particular for what concerns the on-the-fly 
evaluation of scale uncertainties obtained via reweighting 
way as for bottom-associated Higgs production 


56 . This has been done in an analogous 


28 , by splitting the cross section in parts that 


factorise different powers of yb, i. e. yji Vbyt and y^. 

^They appear explicitly separated also in the YUKAWA and MASS blocks of the SLHA cards l55|. 
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3 Results 


In this section, we present four-flavour scheme predictions of charged Higgs boson production 
at NLO matched to parton showers. This calculation has never been performed before in the 
literature. Several differential distributions that are reconstructed from the final state particles in 
tH~b production are studied. We investigate the role of the shower scale in this process, discuss 
the impact of the ybUt interference term and compare our reference predictions at (N)LO-I-PS to 
the f(N)LO results. For matched predictions, both HerwigH — h and PythiaS are employed. We 
conclude this chapter with a comprehensive comparison of 4FS and 5FS distributions, in which the 
effects of higher order corrections, the impact of the choice of the shower scale and the dependence 
of each scheme on the different Monte Carlos are analysed. 


3.1 Settings 

We present results for charged Higgs boson production at the LHC Run H (v^had = ISTeV) by 
considering two scenarios: a lighter (to//± = 200 GeV) and a heavier (to//± = 600 GeV) charged 
Higgs boson. For simplicity, we set tan/3 = 8 throughout this paper. At this value, and yf terms 
are of similar size and the relative contribution of the ytyt term to the total cross section is close 
to its maximum. Results for any other value of tan /3 can be obtained by a trivial overall rescaling 
of the individual contributions according to their Yukawa couplings (//{, by tan/3, y* by 1/tan/3). 
Therefore, we preserve the generality of our results by studying the y^, y^ and ytyt contributions 
separately. 

We show results obtained with the NNPDF2.3 set at NLO and the NNPDF3.0 set at 


LO. To obtain consistent predictions, parton distribution functions (PDFs) computed in the proper 
flavour number scheme are used: we interface our NLO (LO) calculation with the NNPDF2.3 
(NNPDF3.0) with n/ = 4 and n/ = 5 active flavours for the 4FS and 5FS respectively. The 
mismatch between the PDF sets used in the LO and NLO computations is due to the absence of 
a public set of non-QED LO PDFs in the NNPDF2.3 family. This does not affect the accuracy 
of our results, given that the LO PDF sets exhibit a theoretical uncertainty which is larger than 
the difference between the two NNPDF families. The strong coupling constant is consistent with 
as{Mz) = 0.118 for the 5FS NLO parton densities and as{Mz) = 0.1226 for the 4FS NLO ones|^ 
The heavy quark pole masses are set to 


^poie _ ^ |.]^g 4PS), = 172.5 GeV. 


pole _ 


( 8 ) 


At one loop, the value of the bottom pole mass translates into a MS mass 


mbirfib) = 4.3377 GeV. 


(9) 


Finally, our central renormalisation and factorisation scales are set to 

^J■R,F = Ht/3 = ^ X! +PTii)^, 

i 


( 10 ) 


where the index i runs over all final state particles (the top quark, the charged Higgs boson and 
possibly the extra b quark and/or light parton) of the hard process. For vanishing transverse 
momenta of the external particles, our scale choice corresponds to the factorisation scale set in 
the 4FS calculation of Refs. [Toll^. In the following, scale uncertainties are obtained by varying 


fj.F and fj-p independently by a factor of two around their central values, given in Eq. (10). We 
have checked that, particularly for our reference 4FS NLO-I-PS prediction, the dependence of 


^This is the value of ots{Mz) associated with the NNPDF23_nlo_as_0118_nf4 set: the 4FS sets are constructed by 
evolving backwards the 5FS PDFs and the strong coupling constant from the Z mass to the threshold associated 
to the bottom PDF. They are then evolved upwards from the bottom threshold to higher scales by setting ny = 4. 
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the distributions on the shower scale /Xsh, when varied by a factor in the range [l/-\/2, V^], is 
rather mild and significantly smaller than uncertainties associated with the renormalisation and 
factorisation scales; we therefore refrain from including uncertainties associated with jigh in what 
follows. Furthermore, we will not discuss any PDF systematicsj^ 

Jets are reconstructed via the anti-fc^ algorithm [^, as implemented in Fast Jet with 

a distance parameter Ai? = 0.4 and subject to the conditions 

PT(j)>25GeV, |77 (j)I<2.5. (11) 


For fixed-order computations jets are clustered from partonic final states, while in simulations 
matched to parton showers jets are made up of hadrons; b jets are defined to contain at least one 
b quark (at fixed order) or B hadron (in matched simulations). 

In our simulations we keep the charged Higgs boson stable, while we decay the top quark 
leptonically (although the leptons from the decay will not affect any observable we consider) in 
order to keep as much control as possible on the origin of QCD radiation. The task to decay the 
top quark is performed by the parton shower for (N)LO-I-PS runs, while at fixed order we simulate 
the decay t bW in an isotropic way (in the t rest frame) at the analysis levelj^ No simulation of 
the underlying event is performed by the parton shower. 

Let us conclude this section by addressing one further point, which is crucial when processes 
with final-state b quarks are matched to parton showers: the choice of the shower starting scale /ish- 
Such processes are known to prefer much lower values of the renormalisation and factorisation scales 
than the one naively identified as the hard scale of the process (s). In fact, the shower starting scale 
and the factorisation scale emerge both from the same concept, namely the separation of soft and 
hard physics. Furthermore, it has been argued in Ref. 28 for the associated production of a neutral 
Higgs boson with bottom quarks that the shower starting scale (limiting the hardest emission that 
the shower can generate) should be set at similar values, i. e. well below s. Following the arguments 
made in Ref. , we check their validity in the case of charged Higgs boson production. We shall 
stress at this point that the following discussion applies both to our reference scenarios with 
rriH- = 200 GeV and rriH- = 600 GeV, although we refrain from showing explicit results for the 
latter. 

MadGraph5_aMC@NLO assigns a dynamical shower scale chosen from a distribution in the 
rang^ 

0 . 1 . 

-ys < /ish < pS, (12) 


where F is a parameter that drives the bounds of the distribution, and whose default value is 
F = 1. With such a default setting the effective value of /ish, namely the maximum of the /ish 
distribution (which for simplicity we will refer to as just /ish in the following), is indeed much 
larger than /i^.i?- Furthermore, considering the transverse momentum distribution of the Born- 
level “system” (pj,(sys))j^ which is maximally sensitive to the interplay between the fixed-order 
prediction and the shower, the NLO-I-PS distribution (in particular in the 4FS) does not match 
the fixed-order NLO (fNLO) one at large px for F = 1. This can be deduced from Fig. when 
comparing the crosses (NLO-fPS for F = 1) to the solid curves overlayed with points (fNLO). On 
the contrary, we observe a clearly improved high-pr matching of the NLO-I-PS results to the fixed- 

®Note that scale variations due to pp and XR 3-s well as PDF uncertainties are computed at no extra CPU cost 
using the reweighting procedure of Ref. [56| . 

^ Such an approach neglects spin-correlation in the decay of the top quark. However, within the Mad- 
Graph5_aMC@NLO framework, spin correlation can be included in (N)LO-I-PS runs by decaying the top quark 
with MadSpi n [62| . 

®See Ref. 1321 tor further details. 

®Note that the Born-level system is unambiguously defined only in a fixed-order calculation, being in our case 
the charged Higgs accompanied by the final state top and bottom quark. At NLO-I-PS we define it to include the 
hardest B hadron (instead of the bottom quark), which does not originate from the top decay; in this case, MC-truth 
is used. 
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Figure 2: Transverse momentum of the H~bt system for niff- = 200 GeV in the 4FS at fNLO 
(green dotted-solid for the yl term, orange dotted-solid for the term), and at NLO-l-PS with F= 
1 (green dots for the yl term, orange dots for the y1 term) and F= 4 (green solid for the yl term, 
orange solid for the y^ term). We show predictions matched with PythiaS (left) and HerwigH— h 
(right). The insets show the ratio of the curves in the main frame over the fNLO prediction, for 
both the yl and the yl terms. 


order ones by choosing a reduced shower scale corresponding to F = 4 (solid curves)]^ Indeed, 
such a choice brings /igh much closer to the value of the renormalisation and factorisation scales. 
We have also checked that the agreement among PythiaS and HerwigH — h improves (although 
often only marginally) when differential observables in the 4FS are computed with F = 4. 

In conclusion, although for this process we do not reproduce all results of Ref. with the 
same significance, we still find sufficient evidence that F = 4 is favourable in many respects and 
make it our default choice. In Sect. |3.3[ we shall further study the impact of this choice when 
comparing the 4FS and 5FS results: by setting F = 4 an improved agreement between the two 
schemes at the level of shapes is observed. 


3.2 Four-flavour scheme results 


We now turn to our phenomenological results for charged Higgs boson production. Let 
consider state-of-the-art 4FS predictions, which, as will be shown, constitute the most 
differential results for observables exclusive i n the d egrees of freedom of final-state bottom 
We split this section into two parts: in Sect. 3.2.1 we limit our study to the dominant yl 
contributions, while the ybyt contribution is considered in Sect. [X2?^ 


us first 
reliable 
quarks, 
and yl 


9 


























a{mH~ = 200 GeV) [fb] 


NLO 


LO 




2/t 


vl 


yl 



Inclusive 

50 40+^’^ ®^ 

OU.UU_ig g% 

42 43+12.4% 

An 1 9+52.2% 
^^■^^-31.9% 

9Q ^0+36.3% 


F.O. 

45 47+17.5% 

OQ oi +12.2% 

oo 9^+51.9% 

zo.uy_24 g% 

> Ijb 

Pythia8 

43 44+^’7'^^ 

+'J-++_18.4% 

36 67+12 0% 

Of? OI +52.0% 
00.01_3^ 

25 nq+36.1% 

zo.uy_24 7% 


HerwigH—h 

42.64 

36.04 

36.08 

24.61 


F.O. 

1 1 c:c,+10.9% 

q 76+6.5% 
0 -10.0% 

1 1 99+50.4% 
2% 

7 7q+35.0% 
^•'^-24.1% 

> ‘^jb 

Pythia8 

1 9 c c;+15.3% 

iu.oi_i2 1% 

1 1 70+51.2% 
-^-^■''^-31.5% 

Q -I 9+35.6% 

0--L^_24.4% 


HerwigH—h 

11.03 

9.33 

10.09 

7.00 


a{mH 
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0.588+3+3% 
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HerwigH-H- 

0.602 

0.532 

0.498 

0.376 


Table 1: 4FS predictions for total rates (in fb) for tan/3 = 8. 


3.2.1 ul and contributions at NLO+PS 

We begin our analysis by studying total rates for the production of charged Higgs bosons with 
a mass of 200 GeV and 600 GeV in Table We consider three possibilities: the fully inclusive 
case, the case in which we require at least one h jet, and the one in which two or more b jets are 
tagged. All results are given at both LO and NLO accuracy. The cross sections in which one or 
two b jets are required depend on the approximation and Monte Garlo under consideration. We 
thus report separately results obtained at fixed order, with PythiaS and with HerwigH —\-. Any 
quoted uncertainty is due to scale variation, evaluated as detailed in Sect. |3.1[ they are indicated 
only at fixed order and for results matched with PythiaS, since they show little dependence on 
the specific Monte Carlo. Results for and y^ terms are presented separately. Let us summarize 
the conclusions to be drawn from Table [T] as follows: 

• The scale uncertainty of NLO predictions is substantially smaller than that of the LO ones; 
at NLO the scale uncertainty is larger for the y^ than for the yf contribution (~ 15-20% and 
^ 10-15%, respectively), due to the different renormalisation schemes used for the bottom 
and top Yukawa couplings. 

• Because of our default choice of tan /3 = 8, and yf predictions are of similar size at NLO 
(only ^ 15% different); the difference is larger at LO (~ 30%). As a consequence, the K- 
factors are generally different between the yf and y1 terms; for mu- = 200 GeV, the inclusive 
y^ AT-factor is close to 1.2, while for the j/j term the NLO corrections have a larger impact, 

^Our focus here is on the 4FS prediction. However, similar conclusions, if less stringent, can be drawn from the 
corresponding plots in the 5FS. 
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with K « 1.5; for = 600 GeV, NLO corrections are larger in both cases: K « 1.3 and 
if « 1.6 for the and yf terms respectively. The smaller if-factor for the term is a 
consequence of the fact that, by renormalising the bottom Yukawa in the MS scheme, the 
LO predictions already include a class of higher order corrections. 

• PythiaS and HERWIG^—h predictions for cross sections with 6-jet requirements are consistent 
with each other and well within quoted uncertainties; this holds true both at the LO and at 
the NLO. Even the agreement with the fixed-order results is rather good overall and largely 
within uncertainties. Higher 6-jet multiplicities tend to deteriorate the agreement, although 
only to a small extent. 

• The effect of the cuts on the number of 6 jets is quite moderate: given that there is generally 
a hard 6 jet from the top quark decay, the inclusive rate is reduced by only ~ 10 %, when 
at least one 6 jet is required. On the other hand, requiring a second 6 jet reduces the cross 
section by a factor of 4-5. Scale uncertainties slightly decrease at NLO for cross sections 
within cuts. 


We now turn to differential observables in the 4FS. We only discuss results at fixed order and 
matched with PythiaS. Any differences to the matched HerwigH — h predictions will be explicitly 
mentioned below. Besides, the Monte Carlo dependence will be investigated in more detail when 
comparing distributions in the two flavour schemes in Sect. |3.3[ 

The figures throughout this section are organised according to the following pattern: the main 
frame displays the relevant predictions in absolute size as cross section per birj^for the y^ and 
terms. For the sake of readability, they have been separated by reducing the y^ curve by a factor of 
ten. The hxed-order results are shown using boxes, while simulations matched to parton showers 
are displayed as lines: filled green boxes for fNLO, open red boxes for fLO, a black solid line for 
NLO-fPS and a blue dotted curve for LO-fPS. In the first two insets we display the bin-by-bin 
ratio of all y^ histograms in the main frame normalised to the corresponding NLO-fPS result. The 
difference between the two insets are the uncertainty bands, which in the first one are displayed 
for the LO predictions, while in the second for the NLO ones. The third and fourth insets are 
identical to first two, but for the y^ contribution. 

Let us start in Fig. with the transverse momentum spectrum of a 200 GeV Higgs (top left 
panel) and the associated top quark (top right panel). The two plots are in fact quite similar 
and can be discussed simultaneously: the agreement between the hxed-order and the PS-matched 
simulations is close-to-excellent in all cases. Matching to the PS has only minor effects on the shape 
of the distributions, notable only at small transverse momenta. In other words, the resummation 
effects of the shower are extremely small, as can be expected from observables that are practically 
not affected by any large Sudakov logarithm. The NLO corrections mostly affect the normalisation 
and reduce theoretical uncertainties, rehecting the numbers quoted in Tab. [l] We stress here two 
general features regarding the relation between the y^ and y^ terms which shall hold true for all 
subsequent observables under consideration: Hrst, the MS renormalisation makes the uncertainty 
associated with the y^ curves larger, due to the variation of in yj, ~ fhb{pji)] second, t he ra tio 
of the y1 and y^ contributions at NLO is generally flat. More details will be given in Sect. |3.3| 

For a heavier charged Higgs boson with mass mn- = 600 GeV, the lower panels of Fig. 
display a similar behaviour as in the mn- = 200 GeV case. Larger effects due to the PS are visible 
in the Higgs pt spectrum at low transverse momenta. 

We continue our presentation of the 4FS results with the transverse momentum spectra of 
the hardest (pribi)) and second-hardest (pt( 62 )) 6 jet in Fig. see Sect. 3.1 for our jet (and 
6 jet) definition. In this case, we limit our discussion to ttih- = 200 GeV results, since, apart 
from a naturally smaller cross section, the relative behaviour of the different curves is extremely 
similar for a heavier charged scalar. Although the pribi) and PT{b 2 ) distributions in the main 
frame develop a quite different behaviour in terms of their shapes, the ratios in the insets exhibit 


^The sum of all bins is equal to the total rate, possibly within cuts. 
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Figure 3: Transverse momentum distributions of the charged Higgs (left panels) and the top quark 
(right panels) for mu- = 200 GeV (upper panels) and mu- = 600 GeV (lower panels); predictions 
at fixed order (LO red empty boxes, NLO green full boxes) as well as matched to PythiaS (LO 
blue dashes, NLO black solid) are shown in the main frame, for the yl and terms separately. 
The first and second (third and fourth) insets show the ratios of the curves in the main frame over 
the matched NLO prediction, for the y^ (y^) term respectively, together with the LO and NLO 
uncertainty bands. 
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Figure 4: Same as the upper panel of Fig. but for the transverse momentum of hardest (left 
panel) and second-hardest (right panel) b jet. 


comparable features: in all cases, the tail of the spectra is driven by the order relevant to the 
simulation and NLO corrections slightly soften the spectra. In other words, the fixed-order results 
agree rather well with the corresponding PythiaS ones in the tail. On the other hand, close to 
threshold {px = 25GeV), where resummation effects are enhanced, non-showered and showered 
results exhibit sizeable differences, in particular for the hardest b jet. 

Turning now to somewhat related observables in Fig.js]— the transverse momentum distributions 
of the hardest and second-hardest B hadron—, one may expect rather similar features to the 6-jet 
transverse momentum spectra. On the contrary, their pattern is actually very different; the salient 
feature being that showered results are vastly softer than the fixed-order ones and a substantial 
shape distortion due to the matching with parton showers is observed. In fact, even the peak 
of the pt{Bi) distribution is moved by ~ 25GeV towards the left by the shower. However, one 
should bear in mind that we compare bottom quarks at parton level for the f(N)LO predictions 
with B hadrons at (N)LO-I-PS. The observed differences unavoidably lead to the conclusion that 
fragmentation effects become significant for such exclusive observables. Otherwise, the pattern of 
the PythiaS results is very much reminiscent of 6-jet spectra, displaying a slightly harder LO-I-PS 
shape than at NLO-I-PS. Generally speaking, the relative behaviour of the and curves is 
pretty much alike, including the peculiar increase of the f(N)LO cross section towards vanishing 
Pt{B 2 ). Again, we refrain from showing explicit results for a rriH- = 600 GeV charged Higgs boson, 
since the pattern of the various curves turns out to be very similar to the = 200 GeV case; 
the only difference to be pointed out is a slightly reduced gap between showered and fixed-order 
results for mu- = 600 GeV. 

We investigated a vast number of differential observables, the majority of which follows the 
same pattern as illustrated in Figs. and the NLO corrections are rather flat and lie within 
the LO uncertainty bands, shower effects are moderate and become more substantial the more 
exclusive the observable is with respect to the bottom-quark degrees of freedom. Based on our 


13 














































Figure 5: Same as the upper panel of Fig. but for the transverse momentum of hardest (left 
panel) and second-hardest (right panel) B hadrons. 


findings, we do not expect larger effects for 6-jet observables than for distributions relevant to 
B hadrons. Therefore it is worth to consider the invariant mass M{Bi, B 2 ) of and the distance 
AR{Bi, B 2 ) in the rj — (j> plane between the two hardest B hadrons, which are displayed in the 
left and right panels of Fig. [^respectively. The reader should keep in mind that in the fixed-order 
cases the two hardest 6 quarks are used instead. Since there are no salient differences for larger 
Higgs masses we only discuss the mn- = 200 GeV results. 

Although the invariant mass is quite a different observable than the transverse momentum of 
the hardest B hadron, our findings are actually rather similar, but less pronounced: the shower 
substantially affects the distributions, by causing a significant softening whose size exceeds the 
uncertainty bands of the fixed-order calculations. This effect reflects the loss of energy of the 
6 quarks due to their fragmentation into B hadrons. These observations are to a good extend 
independent of the considered contribution (y^ or yf). 

In contrast, the effect of the shower in the AR{Bi, B 2 ) distribution is smaller and becomes 
relevant only at large separations {AR > 4.5). We point out that such large separations are at or 
beyond the coverage edge of the 6 tracking system of LHC experiments, and that if two 6 jets are 
explicitly required these differences are reduced. Effects due to the parton shower are also visible 
at small separations, where secondary g ^ bb splittings can be important. In this case, the two 
B hadrons are likely to be clustered in the same 6 jet. The salient differences between yf and y"^ 
contributions are related to the behaviour of the respective LO curves, which again depends on 
their different relative normalisations; their shapes, on the contrary, are quite similar. 

Let us now conclude our 4FS study by considering jet rates, displayed in Fig. [Tj In the left 
panel, we show jet multiplicities without requiring any 6 tagging, while in the right panel we show 
6 -jet multiplicities. First of all, we recall that in our setup the top quark decays (leptonically) both 
in the shower and at fixed order. For this reason, up to two/three jets and up to two 6 jets can 
appear at fLO/fNLO. Looking at the mu- = 200 GeV results first (upper plots), we can appreciate 


14 
















































Figure 6: Same as the upper panel of Fig. but for the invariant mass of the two hardest B 
hadrons (left panel) and their distance in the rj — (j) plane (right panel). 


the effects of the NLO corrections and of the matching with parton showers. NLO effects have 
the largest impact in the two-jet bin, where the cross section is increased by 35% (65%) for the 
Vb iVt) term. Their effect is minor for lower multiplicities: almost no correction occurs for the 
ul term, while the one is increased by 20%. If in turn we consider the effects of the parton 
shower and compare the fNLO histograms to the NLO-I-PS ones, we infer that the main effect is to 
populate higher multiplicities. Therefore, the zero- and one-jet rates are reduced as a consequence 
of unitarity. Such a reduction is quite important, as the NLO-I-PY8 prediction falls outside the 
fNLO uncertainty band, particularly in the one-jet bin. The two-jet bin is left almost unchanged 
by the shower. We also find that shower effects have a much larger impact at LO, which reflects 
the large uncertainties associated with the LO computations, particularly at fixed order. 

The distribution of events with respect to the number of b jets is displayed in the upper right 
panel of Fig.[^ In this case NLO corrections have the largest effect on the zero- and one-6 jet bins. 
Since the majority of events lies in the one-6 jet bin, the shower moves events from this bin into 
the higher and lower multiplicities, although the effect is in general moderate and within the fNLO 
uncertainties. Only in the zero-6 jet bin, whose rate is however suppressed, the differences between 
fNLO and NLO-I-PS results are larger than in the flavour-unspecific case and the uncertainties 
barely overlap. Overall, the NLO predictions are reasonably close to each other, since higher 
multiplicities (beyond two 6 jets) are phase-space suppressed in the NLO-I-PS simulations. 

For both the jet and 6-jet rates, the fraction of events with jet multiplicities beyond the ones 
already present at the hard-matrix element level (more than three jets and more than two 6 jets) 
is in good agreement between the LO-I-PS and NLO-I-PS predictions, the LO ones being slightly 
enhanced though. At this point, we remark that these multiplicities are utterly Monte Carlo 
dependent and a substantial disagreement between PythiaS and HerwigH — h is found for these 
bins. For the lower multiplicities, however, their agreement is excellent, as will be discussed in 
detail in Sect. 13.31 
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Figure 7: Same as Fig. but for the jet rates (left panels) and &-jet rates (right panels) with 
rriH- = 200 GeV (upper panels) and rriH- = 600 GeV (lower panels). 
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tanS 

Figure 8: Relative contribution of the ytUt term at LO (dashed) and NLO (solid) in the 4FS, with 
respect to the total (N)LO cross section, as a function of tan (3. The two cases mu- = 200 GeV 
(green) and ttih- = 600 GeV (blue) are shown. 


The results of Fig. analysed so far are for rrifj- = 200 GeV. Since these observables are 
particularly relevant for experimental analyses based on jet and 6-jet categories, let us discuss 
explicitly the results for a charged Higgs mass of mn- = 600 GeV, which are displayed in the lower 
plots of Fig. While for jet multiplicities the general features are not so different, some specific 
features are. Apart from a larger Higgs mass changing the LO normalisation, the distribution 
of events at LO appears shifted towards low multiplicities as compared to LO-I-PS, without any 
overlap of the corresponding uncertainties in the first two bins. However, given our findings for 
the lower Higgs mass case, this is expected: the shower shifts events from lower towards higher 
jet multiplicities; this is enhanced for tuh- = 600 GeV due to a generally increased hardness of 
the process. Indeed, the two-jet bin has the largest rate, and the three- and four-jet bins are 
less suppressed than in the lighter Higgs case. NLO corrections slightly improve the agreement 
of showered and fixed-order results, albeit fNLO and NLO-I-PS still fall outside the respective 
uncertainties in the zero- and one-jet bins. 

In the case of 6 jets, on the other hand, the features of the relative curves reflect those discussed 
for the lighter Higgs and no further comments are needed. 

3.2.2 The ytyt contribution 

As we mentioned before, in the 5FS the NLO cross section receives contributions either proportional 
to y^ or to yf. No y^yt term appears, given that it would come from the interference of left-handed 
with right-handed massless bottom quarks. If in turn 6 quarks are massive, as in the 4FS, the ybyt 
term does not vanish any longer, and it is proportional to rn'llQ'^, where Q is some hard scale of 
the process. So far, we have limited our 4FS analysis to the y^ and yf contributions, assuming the 
y^yt one to be suppressed. In this section, we show that this is indeed the case. 

To this purpose, we consider the total cross section for charged Higgs production in the 4FS at 
LO and NLO, and plot the relative contribution —CTj/tyj/cTaii as a function of tan/3 in Fig. with 
CTaii being the sum of all terms. The results are shown for mu- = 200, 600 GeV. The minus sign 
takes into account the fact that the ybyt term is negative. We stress that the ybyt contribution is 
independent of tan /3. As can be inferred from the plots, the relative size of the ybyt term is below 
5% for niH- = 200 GeV, and 0.5% for mn- = 600 GeV. The relative contribution to the cross 
section proportional to ybyt is maximal when the yf and y^ terms are equal, i. e. when 

tan/3^ = ?/(/tan/3^ => tan/3 = 7.27(7.67), (13) 
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Figure 9: Differential comparisons at the between the j/bj/t term and the y^, ones, for ttih- = 
200 GeV. The transverse momentum of the top quark (top left), of the Higgs boson (top right) and 
of the hardest (bottom left) and second hardest (bottom right) B hadron are considered. In the 
main frame the yl (black), y1 (green) and ybyt (blue) distribution are plotted at LO (dashed) and 
NLO (solid), while the inset shows the ratio cr{ybyt)/ (c)?/^) + ^{Vbyt) + ^.t LO and NLO. 
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Figure 10: LO (dashed) and NLO (solid) predictions matched with PythiaS in the 4FS (red for 
yl, violet for j/^) and 5FS (black for light blue for y^), for the transverse momentum of the top 
quark (left) and of the charged Higgs boson (right). Rescaling factors are introduced in the main 
frame for better readability. The first and second insets show the ratio over the NLO prediction in 
the 5FS for the y^ and y^ term respectively, and the scale uncertainty band for the NLO curves. 
The third inset show the differential Ff-factor (NLO/LO) for the four predictions. A charged Higgs 
boson mass ttih- = 200 GeV is considered. 


at LO (NLO), for tuh- = 2OOGeV0 

Let us further investigate the potential impact of the inclusion of the ybVt term on some differ¬ 
ential observables, for such a value of tan /3. In particular, we look at the transverse momentum of 
the Higgs, the top and the two hardest B hadrons for rriH- = 200 GeV, displayed in Fig. From 
these plots we notice that the effect of the y^yt term is peaked at low scales, by reaching at most 
6 — 7% of the full cross section, and is almost the same at LO and NLO. We stress again that 
these numbers have been computed for the value of tan for which the relative ybyt contribution is 
maximal: for larger (smaller) values of tan/3, this contribution is suppressed by a factor l/tan^/3 
(tan^,5) with respect to the y^ (y^) contribution and further reduced for heavier charged Higgs 
bosons. The typical scale uncertainties at NLO (^ 10 — 15%) justify our choice to neglect the 
ybyt contribution in the current analysis. A viable alternative would be to include the relative 
contribution of the ybyt term only at LO, which was shown to be very similar to the NLO one. 

3.3 Four- and Five-flavour scheme comparison 

We turn now to investigate how predictions obtained in the four- and five-flavours schemes compare. 
The two schemes are actually identical up to 6-mass power suppressed terms when computed to 
all orders in perturbation theory, but the way of ordering the perturbative series is different. As a 
consequence, the results in the two schemes may be different at any finite order, while the inclusion 

®The difference between the LO and NLO values is due to the different perturbative order in the running of yf,. 
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Figure 11: Same as Fig. 10 but for the transverse momentum (left) and pseudo-rapidity (right) of 
the hardest 6 jet. 


of higher orders necessarily brings the predictions in the two schemes closer to each other. We 
start by quantifying how the inclusion of NLO corrections improves their mutual agreement. In 
Figs. I0p2 we show, for some relevant observables, the LO and NLO predictions (matched with 
PythiaS) in the two schemes. All figures have the same pattern: a main frame with the absolute 
predictions in the 5FS (black for yl and light blue for y^) and the 4FS (red for y^ and violet for 
y1) at LO (dashed) and NLO (solid). In the first and second insets we show the ratio of the curves 
in the main frame over the 5FS NLO prediction, for the y^ and y^ contributions respectively. For 
the NLO predictions, a band indicating the scale uncertaintjf^ is attached to the curves. In the 


third inset, the four differential NLO/LO AT-factors {yl and yl for 4FS and 5FS) are displayed. 

A general observation is that, as expected, the NLO predictions in the two schemes are much 
closer to each other than the LO ones, in particular as far as shapes are concerned. Differences 
in the overall normalisation reflect the differences in the total cross section, which were already 
discussed in Ref. 22 , while in this comparison we are mostly interested in the shapes. In Fig. 10 


we observe that for the transverse momentum of the top quark and the Higgs boson the difference 
between the two schemes can be compensated by a simple overall rescaling of the total rates 
(cTj^^/cr®^® ~ 0.7) at NLO, while LO predictions in the two schemes have considerably different 
shapes. The same level of agreement should be found also for observables related to the (leptonic) 
decay products of the top quark and the Higgs. Let us recall that in our simulation we do not decay 
the Higgs boson, but we decay leptonically the top quark. The b quark from the top decay mostly 
ends up in the hardest b jet. This explains why the px spectrum of the hardest b jet (left plot 
in Fig. 11) displays a flat ratio between the 4FS and 5FS at NLO, up to ~ 120 GeV. Above that 


value, secondary g ^ bb splittings from hard gluons become more relevant, which is also reflected 
in the growth of the 5FS uncertainty band and IF-factor. A similar behaviour has been observed 


recall that we vary both renormalisation and factorisation scales by a factor of 2 independently about their 
central values. 


20 







































































, 10 '^ 


r.4 

1.2 

1 

0.8 

0.6 

1.4 
1.2 

1 

0.8 

0.6 

2 

1.5 
1 

0.5 


H'bt production at the 13 TeV LHC 

SFSytj^, NLO - 

5FSyt2(x0.1),NLO - : 

LO,NLO+Pythia8, F=4 

4FSyb^, NLO - 

mH-=200 GeV, tanS=8 

4FSyt2(x0.1),NLO - 



5FS yb^. LO (xO.5) . 



5FSy,2 (xO.05), LO . 



4FS yb^. LO (xO.5) . 

- 

. 

4FSyt2 (xO.05), LO . 

- 

anti-kj R=0.4 pj>25GeV'";..'.'; 



ini<2.5 


Ratio over 

5FS, NLO 




NLO/LO K-factor 


20 


40 60 80 100 120 

PT(b 2 ) [GeV] 


H'bt production at the 13 TeV LHC 
LO,NLO+Pythia8, F=4 
mH-=200 GeV, tanp=8 
anli-kj R=0.4 pT>25GeV lril<2.5 


SFSyt,^, NLO • 
5FS yt2(x0.3),NLO 
4FSytj^, NLO 
4FSyt2(x0.3),NLO • 


5FS yb^, LO (xO.5) • 
5FSyt2{x0.15), LO • 
4FS yb^, LO (xO.5) - 
4FSyt^(x0.15), LO • 


1 2 Sfc^Ralioover 5FS. NLO 
1 

0.8 



Figure 12: Same as Fig. 10 but for the transverse momentum (left) and pseudo-rapidity (right) of 
the second-hardest 6 jet. 


in the case of tH production in the SM 29 . The pseudo-rapidity of the hardest b jet (right plot in 
Fig. 111 is mostly dominated by the low-pj^ region, and it therefore also displays a good agreement 
between 4FS and 5FS shapes at NLO. 

Larger differences between the two schemes appear for the second-hardest b jet, which is ex¬ 
pected to be poorly described in the 5FS. In particular, its kinematics in the 5FS at LO is deter¬ 
mined by the shower, while at NLO it is driven by a tree-level matrix element (therefore being 
formally only LO accurate). Predictions for the transverse momentum of the second b jet and 
its pseudo-rapidity are shown in the left and right panels of Fig. [H The 5FS develops large K- 
factors and larger uncertainties, since its LO prediction stems only from the shower evolution. 
Therefore, the 4FS description has to be preferred for these observables, both because of its better 
perturbative behaviour and the proper modelling of the final-state b jets. 

The effects of the different treatment of the bottom quark in the two schemes is even more 
visible for the differential observables related to the hardest B hadron (see Fig. 131. At medium 
and large pt{Bi) and at central r]{Bi) similar effects as for the hardest b jet are observed. At 
variance, the 4FS prediction is suppressed with respect to the 5FS one at low pt{Bi) and at 
large rj{Bi). This is most likely due to mass effects: these kinematical regions correspond to one 
b quark being collinear to the beam. In the 5FS these configurations are enhanced because of 
the collinear singularities, while in the 4FS such a singularities are screened by the 6-quark mass. 
Therefore, even after the PS, the 5FS is reminiscent of the collinear enhancement. In the case of 
the second-hardest B hadron (not shown) these effects are further enhanced. 

Let us make a final remark on the inclusion of the NLO corrections. The NLO/LO AT-factor 
is quite different in the two schemes: in the 4FS the AT—factor appears much more pronounced 
for the ?// than for the term, while in the 5FS it is similar for both contributions. Despite 
that, a remarkable compensation in shape between the LO differential cross sections and the NLO 
corrections takes place, such that the 4FS/5FS ratio at NLO is quite similar for the yf and 
terms. 
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Figure 13: Same as Fig. 10 but for the transverse momentum (left) and pseudo-rapidity (right) of 
the hardest B hadron. 


All the plots discussed so far are relevant to the lighter Higgs under consideration (ttih- = 
200GeV). For a heavier charged Higgs boson (toj^- = 600GeV), the picture does not change 
significantly. The only thing that may be worth mentioning is the fact that, for the second b jet, 
the AT-factor in the 5FS lies much closer to unity than for the lighter Higgs. Such a behaviour may 
be due to the increased weight of the initial-state collinear logarithms resummed by the bottom 
PDFs in the 5FS, which are enhanced at larger masses of the produced particle. Besides, as already 
pointed out in Ref. [23| , collinear logarithms become increasingly relevant the larger the fraction 
of the momentum carried by the initial partons is. 


We continue our analysis by investigating how the choice of the shower scale affects the results 
in the two schemes. We stress once more that our default choice {F = 4) is physically well- 
motivated by the arguments given in Sect. |3.1[ Below we show that this choice also improves 
the mutual agreement between the NLO predictions in the two schemes. A number of differential 
distributions in the 4FS and 5FS for F = 1 and F = A are shown in Figs. mu for both 
rriH- = 200 GeV and 600 GeV. The main frame displays predictions in the 5FS for the term 
(black) and term (red) as well as in the 4FS (in blue and orange respectively). Solid curves are 
used for our reference predictions with F = A, while dashed curves refer to the default choice in 
MadGraph5_aMC@NLO {F — 1 ). The first inset displays the ratio of the curves in the main 
frame over the corresponding ones {y^ or y^) in the 5FS for F = A. The second inset shows the 
ratio of the curves with F = 1 in the 4FS over the corresponding ones in the 5FS. In these two 
insets we can analyse whether F = 1 or F = 4 yields a flatter 4FS/5FS ratio. In the last inset, the 
ratio of the normalised yf and yj distributions is plotted, for the 4FS and 5FS and for F = 1,4. 
The purpose of this inset is to study whether the y^ and y^ contributions develop similar shapes. 

Overall, we observe a smaller dependence on F in the 5FS than in the 4FS distributions. 
A similar behaviour was found in the context of Higgs production in association with bottom 
In any case, the choice of F is almost irrelevant when considering observables that 
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Figure 14: 4FS (blue for the term, orange for the yf one) and 5FS (black for the yl term, 
red for the yf one) distributions for the transverse momentum of the top (left) and charged Higgs 
(right), for m^j- = 200 GeV; the F = 1 (dashed) and = 4 (solid) predictions are shown. The 
first inset shows the ratios of the histogram in the main frame over the 5FS, F = 4 prediction. The 
second inset shows the ratios of the F = 1 histograms over the corresponding 5FS ones. The third 
inset shows the ratio of the normalised y^ histograms in the main frame over the corresponding y^ 
ones. 




Figure 15: Same as Fig. 14 
hardest (right) b jet . 


but for the transverse momentum of the hardest (left) and second- 
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Figure 16: Same as Fig. but for transverse momentum distribution of the top (left) and the 
second-hardest b jet (right) with = 600 GeV. 


are not sensitive to the 6-quark degrees of freedom. In Fig. 14 where rriH- = 200 GeV, we notice 
that charged Higgs and top pt distributions are slightly harder in the 4FS (10-15% at large pr) 
for F = 1. Moreover, the F = A choice yields a remarkably flat 4FS/5FS ratio, much flatter than 
for F = 1. Consequently, F = A improves the agreement between the two schemes. 

This improvement becomes even more visible for observables sensitive to the 6 kinematics. For 
example, Fig. 15 shows the transverse momentum of the hardest (left panel) and second-hardest 
(right panel) 6 jet, where F = 1 provides much harder spectra in the 4FS than F = 4. A similar 
pattern is visible in the case of B hadrons. In this case, however, the general agreement between 
4FS and 5FS significantly deteriorates, as it was already observed. 

In general, similar conclusions can be drawn for rriH- = 600 GeV. The only two distributions 
that exhibit a different behaviour as compared to the rrifj- = 200 GeV case are shown in Fig. [161 
The transverse momentum distribution of the heavy charged Higgs boson (left panel) in the 4FS 
is much more affected by the choice of F, with effects that reach up to 40% in the tail of the 
distribution. In this case the 4FS and 5FS mutual agreement is significantly improved by the 
choice F = 4. The transverse momentum distribution of the second-hardest b jet (right panel) is 
less sensitive to the choice of F than in the lighter-Higgs scenario. 

Jet multiplicities are a class of observables that are particularly sensitive to the excess of 
radiation generated by using F = 1. From Fig. [T^ we conclude that the 4FS generally prefers 
higher jet multiplicities than the 5FS. A similar behaviour was observed also in Ref. 29 and it 


was explained by considering the different colour structure of the initial state in the two schemes, 
besides the fact that the process is generally harder in the 4FS. This tendency is slightly reduced— 
and as a consequence the agreement slightly improves—by the choice of a smaller shower scale 
(F = 4). We remark that the dependence on the shower scale is even less apparent in the heavy 
charged Higgs case, and for 6-jet multiplicities (not shown). 

Finally, the ratio between the and terms (last insets of Figs. Tip?) illustrates that the two 
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Figure 17: Same as Fig. 14 
600GeV (right). 


but for the jet multiplicity with rufj- = 200 GeV (left) and rnn- = 


contributions give remarkably similar shapes for all observables under consideration. Any differ¬ 
ence emerges from the dynamical scale choice at which the bottom Yukawa is computed {rhbifJ-R) 
with fxji = Hx/i, see Sect. 3.1). However, such difference is never larger than 10%, with the yl 
distributions being slightly softer than the ones. 


Finally, we analyse the sensitivity of various observables to the parton shower in the four- and 
five-flavour schemes. To this purpose, we compare results at NLO matched with the PythiaS 
and Herwig-|--I- Monte Carlos. Having verified that the relative behaviour of the two Monte Carlos 
hardly depends on the specific choice of the charged Higgs mass under consideration, we limit the 
discussion to the tuh- = 200 GeV results. 

As already pointed out in the introduction, the parton-shower matching for bottom-quark 
initial states in the 5FS involves some approximations: the initial-state backward evolution is 
based on leading-log accurate gluon splittings and requires the reshuffling of massless into massive 
bottom quarks. For these reasons, the 5FS predictions are extremely sensitive to the specific 
treatment of bottom quarks in a given Monte Carlo. This is most remarkable for 6 -jet/i?-hadron 
related observables, such as the transverse momentum distribution of the second-hardest B hadron 


(displayed in the left panel of Fig. 18). In the 5FS, the HerwigH — h prediction displays a significant 
shape distortion, both at small and large values of the transverse momentum, falling outside the 
PythiaS uncertainty bands at small px- A similar discrepancy is evident in the forward region of 
the rapidity spectrum, shown in the right panel of Fig. |18[ In this plot, the 5FS prediction is larger 
for HerwigH — h than for PythiaS. For both observables, the 4FS results display a significantly 
smaller Monte Carlo dependence. 

An even more spectacular example is provided by the rj — 4> distance between the two hardest 
B hadrons. The results are shown in the left panel of Fig. The HerwigH — h prediction in 
the 5FS features a much higher tail at large separations. This can be traced back to the fact 
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Figure 18: Transverse momentum (left) and pseudo-rapidity (right) distribution of the second- 
hardest B hadron. NLO curves matched with PythiaS (solid) and HerwigH —h (dashed) are 
shown, in the 5FS (blue for the yl term, dark-red for the yf one) and 4FS (green for the y^ term, 
light-orange for the y^ one). The two insets show the ratios of the histogram in the main frame over 
the 5FS prediction matched with PythiaS, for the y1 and y"^ terms separately. Scale uncertainty 
bands are shown around the PythiaS predictions. 


that HerwigH —h tends to produce B hadrons much closer to the beam line when simulations are 
performed in the 5FS. The same behaviour was observed in Ref. 28 , in which it was pointed out 


that such effects are however not relevant when realistic cuts on the B hadrons are imposed, for 
example, when h jets are required. As a matter of fact, we observe a neatly improved agreement 
among all curves in Fig. [I^ when requiring at least two b jets (right panel). 

Looking further at jet (left panel) and 6-jet (right panel) multiplicities in Fig. 20 we observe 
instead that the two flavour-schemes display a similar Monte Carlo dependence. In the flavour- 
unspecific case this dependence is quite small for all jet-multiplicity bins. On the other hand, for 
6 jets HerwigH —I- and PythiaS are in good agreement up to one 6 jet in the 5FS and two 6 jets 
in the 4FS, that is up to the multiplicities described by the hard matrix element at NLO. The 
two 6-jet bin in the 5FS—described only at LO by the hard matrix element—is affected by larger 
discrepancies between the two Monte Carlos, which exceed the uncertainty bands. Larger 6-jet 
multiplicities, generated only at the shower level via gluon splittings, are affected by discrepancies 
as large as 100%: HerwigH —h predicts less 6 jets than PythiaS. 

Generally speaking, the difference between the Monte Carlos is smaller for the 4FS than for 
5FS predictions. Indeed, the 4FS has more differential information at the matrix-element level, 
which reduces the effects of the shower. The only case worth mentioning, in which the Monte Carlo 
dependence is larger in the 4FS than in the 5FS is the distance between two hardest 6 jets at large 
separations. This observable is closely related to the distance of the two hardest B hadrons when 
two 6 jets are required (see right panel in Fig. 19). For both observables, at large separations, the 
4FS predictions matched with HerwigH —h lie very close to the 5FS NLOH-PS predictions, while 
the matching with PythiaS yields a higher tail. However, we point out that such discrepancy 
between HerwigH —h and PythiaS in the 4FS barely exceeds the scale uncertainty bands. 


26 




































AR(Bi,B2) AR(Bi,B2) 


Figure 19: Same as Fig. but for the r] — cj) distance of the two hardest B hadrons with no cuts 
(left) and requiring at least two 6 jets (right). 
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Figure 20: Same as Fig. 18 but for the jet (left) and 5-jet 


(right) multiplicity. 
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Finally, we remark that similar conclusions can be drawn for the heavier charged Higgs case for 
most of the studied observables. However, a heavier charged Higgs reduces the differences between 
the two Monte Carlos (in particular in the 5FS) in the transverse momentum distribution of the 
two hardest B hadrons/& jets. 


4 Conclusions 


We have presented predictions for the production of a heavy charged Higgs boson in a type-H 
2HDM, by explicitly considering a lighter Higgs scenario {rriH- = 200 GeV) and a heavier one 
{rriH- = 600 GeV). Our predictions have been presented for tan/3 = 8, but they are applicable to 
any tan /3 value through a simple rescaling. Furthermore, these results can be straightforwardly 
extended to a type-I 2HDM by rescaling the Yukawa couplings. Details are given in Sect. 6 of 
Ref. 
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For the first time, a fully differential computation has been performed in the 4FS at fNLO and 
NLO+PS accuracy. We have exploited the automatised MadGraph5_aMC@NLO framework, 
by obtaining the NLO version of the 2HDM via the NLOCT package. The model has been 
supplemented by the computation of the bottom Yukawa coupling in the MS scheme, which has 
the advantage with respect to the default on-shell scheme of resumming large logarithms of rrib/^ji. 

Our results indicate that a reduced shower scale with respect to the default one in Mad- 
Graph5_aMC@NLO improves the matching between parton shower and fixed-order results at 
large transverse momenta. For this scale choice, we discussed the effects of incorporating NLO 
corrections and the matching to the parton shower in the 4FS simulations by considering a number 
of differential observables. We found NLO corrections to be generally flat with our choice of the 
renormalisation and factorisation scales. Effects due to the parton shower are important in the 
Sudakov-dominated regions (jets and B hadrons at low px), or when observables sensitive to the 
6-quark fragmentation are considered. These observations have been made separately for the y'^ 
and terms. On the other hand, we argued that the ytyt contribution, appearing only in the 
4FS, can be safely neglected, since its size is smaller than ^ 5% of the total cross section (for 
tan /3 = 8 and rriH- = 200 GeV) and is well within the scale uncertainty of the computation. For 
larger Higgs masses or different values of tan /3 it is further suppressed. 

Besides discussing the new results in the 4FS, we provided a comprehensive comparison to 
the ones in the 5FS, consistently generated within MadGraph5_aMC@NLO. The inclusion of 
NLO(-I-PS) corrections in the two schemes improves their mutual agreement at the level of shapes. 
This agreement follows, although to a minor extent, from the reduced shower scale choice. Differ¬ 
ences remain, however, and they are particularly sizeable for observables related to b jets and B 
hadrons. Given these differences, it was vital to carry out an unbiased analysis of our results in 
order to acquire the most reliable predictions for this class of observables. The proper simulation 
of the signal will be crucial for the experiments to fully exploit the potential of the data collected 
in charged Higgs searches at the LHG. 

Our final recommendation is to use 4FS predictions for any realistic signal simulation in exper¬ 
imental searches. This recommendation is backed by two sets of evidences: first we have proven 
that, for a large number of observables, the 4FS prediction provides a better description of the 
final state kinematics; second, it reduces the systematic error related to the usage of a given par- 
ton shower. Moreover, when matching the NLO calculation to the shower, we recommend to use 
a lower shower scale (by setting F = 4 in our case). This choice provides a better matching 
to the fixed-order computation at large transverse momenta, slightly reduces the parton shower 
dependence and improves the agreement of four- and five-flavour scheme computations. 

Any user interested in the simulation of charged Higgs production with MadGraph5_aMC@NLO 
is strongly encouraged to contact the authors. 
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